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^ • Abstract 

o 

o . 

The interplay of topological constraints and persistence length of ring 

^ i polymers in their own melt is investigated by means of dynamical Monte 

•/^ i Carlo simulations of a three dimensional lattice model. We ask if the re- 

^*~ \ suits are consistent with an asymptotically regime where the rings behave 

C^ . 

Q^ , like (compact) lattice animals in a self-consistent network of topological 

2 ! constraints imposed by neighbouring rings. Tuning the persistence length 

,_u , provides an efficient route to increase the ring overlap required for this 

Q , mean-field picture to hold: The effective Flory exponent for the ring size 

O 

decreases down to z^ ~ 1/3 with increasing persistence length. Evidence 

rS ' is provided for the emergence of one additional characteristic length scale 

dt oc A^ , only weakly dependent on the persistence length and much larger 

than the excluded volume screening length ^. At distances larger than dt 
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the conformational properties of the rings are governed by the topological 
interactions, at smaller distances rings and their linear chain counterparts 
become similar. (At distances smaller than ^ both architectures are iden- 
tical.) However, the crossover between both limits is intricate and broad 
as a detailed discussion of the local fractal dimension (e.g., obtained from 
the static structure factor) reveals. This is due to various crossover effects 
which we are unable to separate even for the largest ring size (N = 1024) 
presented here. The increased topological interactions also influence the 
dynamical properties. Mean-square displacements and their distributions 
depend crucially on the ring overlap and show evidence for the existence 
of additional size and time scales. The diffusion constant of the rings goes 
down from effectively D^ cc A^^^-^^ for flexible rings with low overlap to 
Dj\j oc N~^'^^ for strongly overlapping semi-flexible rings. 

PACS numbers: 61.25.Hq, 61.41. +e, 83.10.Nn, 83.20.Fk 



I. INTRODUCTION 

Unconcatenated and unknotted rings in their melt are relatively compact [|l| . This was 
found in recent computational studies [^|-|^ and was expected on theoretical grounds [^j-Q. 
Qualitatively the squeezing of the rings was attributed to the topological constraints (see 
Fig. ^ (a)) [|T^]. This is in line with much older observations showing that dilute rings repel 
each other much more strongly than their linear chain counterparts due to the entropy 
loss associated with the unconcatenation constraint preventing the two rings to thread 
each other [^. While the usual excluded volume interaction is screened out at high 
chain overlap (i.e., if the ring size R is larger than the size ^ of the excluded volume blob) 



I^,[T3|, topological interactions are expected to dominate the conformational properties 
if the number of overlapping rings p ~ R^(f)/N [N being the chain mass, the monomer 
density) becomes high enough. 

In a simple Flory-like argument Gates and Deutsch (CD) 0] argued that the number 
of degrees of freedom lost for a typical ring due to its topological interactions with 
neighbouring rings increases like p°' where a is an unknown exponent. This free energy 
term, which favors decreasing the ring size, has to be balanced by the entropy penalty 
for squashing a ring. To be specific, this penalty was assumed to be the same as for 
a Gaussian chain. Adding these contributions and minimising over R oc N" yields the 
Flory exponent u = {a + l)/(3« + 2). The simplest possible estimate for a is to say that 
that roughly one degree of freedom is lost for each of the p neighbours which the ring is 
prevented from threading, i.e. a = 1 and hence u = 2/5. This is very close to the value 
u ~ 0.39 found by us in our previous study IQ. Note that the overlap number p oc N^^^ 
increases extremely weakly and that the only intrinsic length scale in the CD picture is 
the ring size R itself (possibly renormalised in terms of excluded volume blobs of size ^). 
Again this was found to be qualitatively in good agreement with the simulations reported 
in our first paper |Q. From the CD picture one expects a similar density crossover scaling 
for rings as for linear chains: The ring size R{4>) reduced by the ring size Rq of the dilute 



reference ring should scale with 0/0* where 0* ~ N/Rq is the crossover density. Having 
focused in ref. Q on flexible rings at one fixed density we were unable to verify this 
implicit scaling assumption of the CD picture. 

A different picture comes from the extensive studies on (isolated) ring polymers in gels 
(see Fig. |T] (b)), often modelled by so-called lattice animals (LA), depicted in Fig. |T] (c) 
PJ^J^. If the ring size R becomes larger than the typical distance dt between the fixed 
topological obstacles (represented as squares in the figure) the rings are forced to retrace 
their paths and the fractal dimension df = l/u becomes that of a strongly branched 
object. The question is now if it is possible to use this well-understood model via a 
standard mean-field argument also for strongly overlapping systems of rings in their own 
melt. Are the surrounding rings able to generate (in a self-consistent manner) around 
a reference ring a fixed mesh of topological obstacles? And if so, is dt (x R — then LA 
and CD picture would be indistinguishable from the scaling point of view — or does it 
introduce an additional ring length independent scale? In the latter case, obviously, the 
above mentioned density crossover scaling would not work. (However, one can strictly 
only expect this LA picture to hold in the high entanglement limit.) In any case, there 
is a catch: While increasing p oc b^N^"'^, b being the persistence length, the effective 
exponent z/(p) is supposed to drop down to 1/3 (and even down to z/ = 1/4 for an ideal 
Gaussian LA within an intermediate range of ring sizes A^) . Hence, the LA picture 
might not be self-consistent as already stressed in ref. 0. It is crucially the prefactor 
b^ (and not the inefficient N^'^^p^~^) which allows the simulator to control the overlap. 
(This assumes that the persistence length only weakly affects the hypothetical length dt 
which has to be checked a posteriori.) . 

This is the route we have taken in this study to test the scaling predictions (rather 
than the variation in density which we plan to investigate in a subsequent study ||14|| ). 
Indeed, the persistence length b turns out to be a very efficient way to vary the overlap 
number p (additional computational overhead and reduced diffusion constants taken into 



account) compared to the ring mass (which we have however increased from A^ = 512 to 
A^ = 1024). This allows us to more severely put a test on the scaling predictions of the 
CD scenario than we were able in our previous study 0]. In contrast to that work, the 
new evidence presented here shows the emergence of (at least) one further length scale 
alongside R which we identify with dt- Tentatively, our data are consistent with a broad 
crossover towards the LA picture in the limit of high chain overlap which is attained 
by increasing b so that the chains become more extended (though more compact in the 
scaling sense of smaller z/) . A typical semi-flexible coil at higher overlap is presented in 
Fig. ^ on the right. 

Of course, the CD and LA pictures are not necessarily contradictory in that they might 
provide useful heuristic descriptions in different overlap limits. In the computationally 
important regime (depending on b) where the rings expel successfully neighbouring chains 
and form relatively dense coils, the mean-field assumption of the LA approach has to 
break down. A typical chain for low overlap is shown Fig. on the left. The CD picture 
is a priori a good candidate to describe this regime where overlap p and topological 
interactions are weak. A possible choice for the unknown CD exponent in this regime 
is the limiting case a ^ and, hence, u ~ 1/2, i.e. the topological constraints mainly 
contribute logarithmic corrections to a closed |jlO| Gaussian chains of blobs. Evidence for 
this (justifying a posteriori the squeezing term used by CD) is presented below. 

Our paper is arranged as follows: In the next section |I| we give a short synopsis of 
the model and simulation technique used. To understand the special effects linked to 
the unconcatenation constraint we need reference date at the same chain and persistence 
length to compare with. This is provided in section |T| where we review some properties 
of dilute semi-flexible ring polymers. In the subsequent section ^ we investigate the 
statistics of ring polymers in their melt as a function of the stiffness, comparing them 
with dilute rings and dense linear chains from ref. |T5|,|T6|. (Note that in most figures 



we compare features discussed subsequently in section |T| and |IV|.) Possible scaling 



scenarios are discussed and a detailed analysis of the local fractal dimension is presented. 
The following section |V| presents briefly our first results on the dynamics of semi-flexible 
ring polymers. We conclude with a summary of our results. 

II. ALGORITHM AND PARAMETERS 

As in our previous study we investigate the properties of unknotted and unconcate- 
nated rings in the framework of the bond fluctuation model (BFM) [|l^]. Many static and 
dynamic properties of linear chains are known for this computationally efficient, coarse 
grained lattice model. A small number of chemical repeat units is mapped onto a lat- 
tice monomer such that the relevant characteristics of polymers - excluded volume and 
connectivity - are retained. Each monomer blocks a unit cell of the three dimensional 
cubic lattice from further occupancy. Adjacent monomers along a polymer are connected 
via one of 108 allowed bond vectors of lengths 2, ^/5, y/6, 3 and vTO. Here and in the 
following all spatial distances are measured in units of the lattice spacing. This set of 
bonds allows for 87 different bond angles and, hence, results in a good approximation to 
chains in continuous space. The bond vectors are chosen such that the local excluded 
volume interactions prevent the rings from crossing each other during their motion. This 
conservation of the topology ensures that the rings remain neither knotted with them- 
selves nor concatenated with one another during the course of the simulation. We evolve 
the ring conformations via random local monomer displacements. 

As explained in the Introduction, a crucial parameter is the overlap number p, which 
increases only very weakly with A^, but much more strongly with the persistence length. 
In order to tune the persistence length, we impose a simple intramolecular potential which 
favours straight bond angles: E{Q) = a cos(G). Here denotes the complementary angle 
between two successive bonds and a the dimensionless energy scale (setting fc^T = 1). 
This potential has been used and investigated in various studies on linear polymer chains 
at ffiling fraction = 0.5 of occupied lattice sites [|IB],|rU|. This is of importance because it 
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is reasonable to assume that on short length scales (for long enough chains) neither ring 
closure nor topological constraints are pertinent. This statement will be corroborated 
later in section |V| when we discuss the local ring structure (Fig. ^). Note for now that 
other local quantities like the mean bond length /, (cos(B)) (i.e. the mean stiffness 
energy per monomer) or the acceptance rate A are identical for rings and their linear 
counterparts. Hence, one expects to find the same local rigidity b{a) for rings as for linear 
chains where 6 = 6^ is easily obtained from the mean-square end-to-end vector {RD and 
the known Flory exponent of linear chains. See Tab. | and Tab. |l| for the dilute (0 = 0) 
and the dense limit (0 = 0.5) respectively. The systems containing dilute rings and linear 
chains have been simulated for reasons of comparison. Not surprisingly the quantities 
featured in Tab. | and Tab. || depend somewhat on the volume fraction [|15|,|16[. Note 



that at the density = 0.5 of occupied lattice sites, many static and dynamic features 
of molten linear polymer materials are reproduced by the BFM. For example, the single 
chain conformations obey Gaussian statistics down to the screening length ,^(o" = 0) ~ 6 
of the excluded volume interaction obtained by the static structure factor |T^. 

We have restricted ourselves to values a < 3 because we found tentative evidence 
for the onset of a nematic order at values cr > 4. Due to the lower conformational 
entropy rings tend to order nematically at lower stiffnesses than linear molecules. (For 
linear chains an isotropic to nematic transition is found close to a ~ 6 [0.) Moreover, 
increasing the stiffness also decrease the number of statistically independent segments, 
and we want to keep finite A^ effects low, even for small ring lengths. 

The algorithm described above has been implemented on a massively parallel CRAY 



T3E computer |jT8[. Using a 2 dimensional geometric decomposition of the simulation 
grid of linear extension L = 128, we employ 64 T3E processors. The simulations involve 
about 40 000 hours of single processor CPU time. 

As shown in Tab. fll|, this allows us to investigate ring systems with 131072 monomers 
and chain masses up to A^ = 1024. The data for flexible rings have been used as starting 



configurations for simulations at liiglier persistence lengtli. In most cases rings diffused 
at least a spatial distance of their radius of gyration before any conformational data were 
collected. The simulation runs were extended up to five times the relaxation time for 
sampling the ring conformations. Note that we could not meet this stringent criterion 
for the semi-flexible {a > 0) systems at our largest mass A^ = 1024. The resulting data 
appear to be time independent, and we presume them to be equilibrated, but they have 
not diffused over a radius of gyration. These three data points (a = 1, 2 and 3) have to 
be taken with care. 

As can be deduced from Tab. |T| a relatively small increase in persistence length is 
a very effective way to increase the overlap number p oc Rj' /N . To give some numbers, 
increasing the stiffness parameter a from to 3 amounts in increasing p by a factor 3 
for A^ = 256 (seen in Fig. Q) and by a factor 5 for A^ = 512. To achieve a similar effect 
by tuning the contour length requires an increase by two, respectively three, orders of 



magnitude |^| 



III. DILUTE RINGS REVISITED 

As a reference for the subsequent section on strongly overlapping rings (where the 
unconcatenation constraint matters) we briefly revisit some properties of flexible and 
semi-flexible dilute rings and their linear counterparts. Basically, the remaining unknot- 
tedness constraint and the requirement of ring closure [|10| do not matter: Dilute rings 
behave broadly like linear chains. 

The simplest quantity to characterise the chain structure as a function of the stiffness 
parameter a is the overall chain size versus chain length A^. The size of the rings is 
measured first with the usual mean-square radius of gyration (R'i). As a second measure, 
we sample the mean-square distance (-R^) = {{Ri — Ri+nY) (averaged over all monomers 
i) between pairs of monomers that are n monomers apart along the contour. In particular 
we define for rings the mean-square ring diameter {R^) = {Rn)\n=N/2- (For the linear 
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chains (i?g) = {E^\n=N denotes the usual mean-square end-to-end distance.) In Fig. ^ 
the diameter R^. = {Rl) for flexible and semi-flexible dilute rings is plotted versus A^. 
Not included are the linear counterparts of same A^ and a. (But see Fig. ^) The data 
are well fitted for all the persistence lengths considered by Re = bgN'^ with the classical 
Flory exponent u = uq ^ 3/5 |2T|. The effective bond length be encodes the persistence 
length effect. It is tabulated for linear chains ("L") and rings ("R") in table |l[ Similar fits 
have also be performed for the radius of gyration. We find, independent of chain size and 
persistence length, the following relations: 6^/&f ~ 1.89, bg/b^ ~ VQ and bf/bf ^ 1.79. 
Note that in fitting we have disregarded the three smallest masses A^ = 16, 32 and 64 to 
minimize finite size effects. We have checked for (the remaining) finite size effects (not 
shown): The (negative) curvature in the data points is surprisingly weak and similar 
for linear chains and rings. It is slightly stronger for the radius of gyration than for the 
diameter which probes larger distances. While we cannot rule out that the asymptotic 
exponents of the ring chains are slightly different from their linear counterparts, for all 
practical purposes of relevance here we may conclude that neither the ring-closure nor 
the non-selfknottedness constraint are pertinent. 

We are now in the position to rescale the effects of the persistence length for both 
linear chains and rings by estimating for the linear chain the number gk of correlated 
monomers along the chain and the associated Kuhn length k = Ig^ {I being the mea- 
sured mean length of the BFM bond). Rewriting R^ = beN" = k{N/gkY yields 
dk = (be/iy^^^"'^^ and k = l{be/lY^^^~'^^ ■ Obviously this scaling-based definition of k 
and gk is arbitrary to within prefactors of order one. The proposed rescaling works suc- 
cessfully for dilute linear chains and rings as shown in Fig. ^. Here we have plotted the 
reduced diameters u = Re/k versus the number of statistical segments v = N/g^. A 
similar figure was obtained for the radius of gyration (not shown). Note that only a 
part of the linear chain configurations used to characterise the stiffness effects for dilute 
systems are included in the figure. 



In Fig. ^ we show the bond-bond correlation function < ejei+„ >, where e, denotes 
one of the bond vectors of monomer i (say the "right" one in a given hst). Only values for 
rings are included (the values for linear chains being identical for long enough chains). As 
it should, the bond-bond correlation functions for different persistence lengths collapse 
onto a single scaling curve, when plotted versus n/gk (not shown). Note that this scaling 
and the (supposed) exponential decay is sometimes preferred to define the number of 
correlated monomers Qk and the persistence length k rather than our definition based on 
the measured effective bond length b and the well-known (and for rings confirmed) Flory 
exponent. However, the correlation function for dilute chains (both linear and rings) 
does certainly not decrease as a pure exponential (not shown). This is due to (non- 
universal) short range packing effects and (more importantly) to long range excluded 
volume correlations. In addition to this we prefer our method for statistical reasons. 

Additional information about the local structure of the flexible and semi-flexible rings 
comes from the structure factor S{q) which is of direct experimental relevance. The slopes 
of the structure factor in log — log coordinates define the differential fractal dimension 
dfil) = —d^og{S{q))/d\og{q). The general shape of S{q) was discussed in ref. [^. We 
present here only df{q) which contains all the information. It is equivalent in the limit 
of long chains {q -^ 0, qRg ^ 1) to the inverse of the Flory exponent v. Fig. ^(a) shows 
df{q) for dilute rings (thin symbols) of mass A^ = 1024 plotted versus 5r = 271 /q. For 
large distances the rings are well described by the classical Guinier expansion S{q) ~ 
A^(l — {RgqY /?>) (not shown); hence df{qRg — > 0) ^ 0. For small distances 6r ^ / one 
finds the usual spurious Bragg-peak of the structure factor [Q due to local packing effects 
giving rise to negative (^/-values. We note a weak cr-dependence due BFM related lattice 
effects: A slightly different set of bond vectors is preferred. This is in line with the fact 
that also the mean bond vector varies weakly with persistence length as shown in Tab. |. 
In between both limits we find a broad plateau for the flexible rings with df = 1/z^o ~ 1.7 
(upper horizontal line) due to the excluded volume interactions. Semi-fiexible rings show. 
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as expected from linear chains, a shoulder at about df = 1 (upper horizontal line) due to 
the local rigidity of the rings. Note that the curves for linear dilute chains (not included) 
are very similar. The only qualitative difference is the "hump" between the plateau 
and the Guinier regime. It is due partially to the ring closure as one can easily see by 
calculating the static structure factor for a Gaussian ring (solid line). 

IV. CONFORMATIONAL PROPERTIES OF DENSE RINGS 

A. Compact, but strongly overlapping rings 

While dilute rings essentially behave like their linear counterparts, this becomes very 
different for unconcatenated rings in the melt. This is shown in the dramatic decrease 
of the chain size of rings (full symbols) in Fig. |^ as compared with dilute reference rings 
(open symbols). Certainly, linear chain sizes also decrease with increasing density due 



to the screening of the excluded volume correlations |12,15|, but to nowhere near this 
extent, as Fig. ^illustrates. Moreover, the slopes of i?e versus v = N/g^ decrease strongly 
with persistence length. (See Tab. |T[) For flexible rings {a = 0) the data points (for 
N = 128 to A^ = 1024) are fitted by an (effective) exponent u = 0.4. For our stiffest 
rings we observe z/((j = 3) ~ 1/3. This last result is indeed consistent with the compact 
LA scaling behaviour predicted for rings at high enough overlap p. 

These results are also in agreement with the density distribution of a ring around its 
center of mass, shown for rings in Fig. ^ We see that for flexible rings (main figure) the 
overlap p must be very small; nearly all monomers from neighbouring rings are expelled. 
From the Flory exponent u ^ 0.4 measured above one expects the density inside the ring 
to decrease as p ~ A^~°-^. This is confirmed. In the inset we present the corresponding 
plot for semi-flexible rings {a = 3). Since the effective Flory exponent is z/ = 1/3 we 
expect the density inside the ring to be independent of the ring size. The ring length 



independence of the density is observed for A^ ~ 128 for a = 2 (not shown) and A^ ~ 64 
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for (7 = 3. This supports the idea that the LA regime is reached earher for stiffer rings. 
Though rings in this regime show compact scahng (z/ ^ 1/3) the asymptotically constant 
density p is fairly small. Thus the density of oi/ier rings in the correlation hole of a given 
ring is larger ioT stiff than for flexible rings (Fig. |^). This holds for the ring sizes studied; 
it does not exclude similar strong overlap behaviour for flexible, but much larger rings. 

B. First evidences for additional length scale 

We have to stress, however, that the above exponents are only effective exponents. 
The curves are strongly (negatively) curved as one can easily visualise by plotting R/N^^^ 
versus A^ (not shown). Alternatively, one can characterise this curvature via the differ- 
ential Flory exponent i^(A^) defined from the increase between chain length A^ and 2N 
1^. This reveals (not shown) that z/(A^) decreases continuously from u ^ 1/2, i.e. nearly 
Gaussian behaviour, for smaller rings down to the slopes indicated in Fig. |^ for the largest 
masses we were able to simulate. In our previous study ^ we have attributed the ob- 
served curvature to classical finite-size effects which also appear for linear chains. That 
is, we have not attributed them to topological effects, but excluded volume interactions 
visible due to the finite number of blobs. If one admits this as the only physical origin of 
the curvature (excluding, e.g., additional length scales) it is indeed reasonable to attempt 
to obtain the asymptotic exponent by classical finite-size extrapolation method. In this 
method, the local slope i^iN) (obtained from the diameter and the radius of gyration) 
is plotted versus the reduced blob size ^/R where R is itself given in a self-consistent 
manner by the asymptotic behaviour. This works generally well for linear polymers. 
Proceeding along these lines we obtained an asymptotic exponent oi u ^ 0.4 (which we 
now obtain directly after including A^ = 1024). This value happened to coincide with 
the CD picture (a = 1 |0) which made us perhaps more confident than we should have 
been (in ref. [^) of having characterised the asymptotic behaviour. 

The excluded volume effects (and their screening) do certainly contribute to the gen- 
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eral crossover scenario. This is shown below when we discuss the structure factor (Fig. P). 
But a danger remains that we may have missed additional physics, in the form of an ad- 
ditional length scale (not included in the simple CD picture), which controls the behavior 
of molten rings. 

Thus one of the key claims made in ref. [^ was that the only relevant length scale in 
dense rings is the chain size itself. This was shown in various scaling plots, e.g., of the 
reduced mean ring "area" (|a|)/(i?p = const. The area a is defined as a signed quantity 
(the component in that direction of the vector area of a ring) that vanishes for any 
configuration in which the ring exactly retraces its own steps. Qualitatively, one expects 
that a ring which is only very weakly threaded by the surrounding rings has a very small 
area. This should hold in particular while we approach the LA limit. The reduced area 
(|a|)/(-R^) has been (log-linear) plotted in Fig. ^ versus ring mass A^ for dilute rings 
(empty symbols) and rings in the melt (full symbols) of different flexibility. For (long 
enough) dilute rings the ratio is indeed chain length (and persistence length) independent. 
For the dense limit however, we now find evidence for (logarithmic) corrections to the 
expected CD scaling (plateau). The reduced area decreases systematically as a function 
of mass — this is in contrast to what was deduced from a (linear-linear) plot in ref. [Q. 
Note that the ratio decreases also with respect to a: stiffer rings are likely to get less 
threaded by other rings. This is consistent with the LA picture where stiffer (and hence 
larger) rings are forced to retrace their own steps in a network of fixed obstacles. 

C. Failure of classical scaling analysis without additional length scale 

More stringent tests for the CD picture are posed by the following two scaling analyses 
which try to allow for a varying persistence length under the assumption of no new 
additional length scale for topological interactions. They are presented in Fig. ^ and in 
the inset of Fig. ||. 

In the first we replot the reduced ring size u = Re/k versus the number of statistical 
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units V = N/gk as we did in Section |n| for dilute rings. The values of k and gk used in 
the Fig. ^ for the dense linear chains and rings are the same as for the dilute systems ||22|| . 
The difference to the behaviour of linear chains is striking: The ring data for different a- 
lines diverges while the linear chain data clearly collapses (on the expected slope z/ = 0.5). 
This is a physically unsound result; flexibility- dependent universality classes for rings are 
difficult to accept. 

The second scaling test is the "classic" u = R{a, = 0.5)/-Ro versus v = 0/0* ~ N/R^ 
mentioned in the Introduction [|12|. This is shown in the inset of Fig. p. Both ring 



diameter (full symbols) and radius of gyration (empty symbols) at = 0.5 have been 
included. As reference chain size Rq we used the measured radius of gyration for a dilute 
ring of given a. Again, the scaling attempt fails to account for data of more than one 
chain stiffness. 

For the moment we may conclude that there is clear evidence for one or more addi- 
tional length scales. At this point, we still do not know how this hypothetical length scale 
df depends on ring and persistence lengths, nor do we know what physics it represents: 
however, since it arises for rings and not for linear chains, we can presume a topological 
origin. 

D. Evidence for a chain length independent length dt 

Another striking effect which we are able to see due to the persistence length variation 
is shown in Fig. ^ for the bond-bond correlation function of dense rings. As before for 
dilute chains this correlation function becomes chain length independent for the large 
masses indicated in the figure. For fiexible rings the correlation function drops down 
slightly below zero within three monomers and approaches then zero from below. When 
the persistence length is increased a most remarkable negative correlation becomes visible 
indicating that the polymer is likely to fold back after 10 monomers for {a = 3) | |23|| . 
Needless to say that the correlation function for dilute rings does not show anything 
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faintly similar. The position and the depth of this anti-correlation dip increases with a. 
We stress that the position and the shape of the dip are chain length independent. If 
this effect has something to do with the additional length scale dt mentioned above, as 
we believe, this is a piece of evidence for its chain length independence. 

E. A new scaling scenario 

This observation forms the basis for the next scaling proposal attempted in Fig. ^ 
We assume here that there is one additional, mass independent length scale dt- Hence, 
the associated number gt of monomers between the topological obstacles is also mass 
independent, but depends on the local conformational properties. These are complicated 
as will be revealed below (see Fig. ^ and Fig. |10D. However, as we have observed above 
(Fig. ^, short rings are reasonable approximated by Gaussian statistics and we may 
write Qt = Qk (dt/k) ''^^ where ui is some effective exponent close to 1/2 characterising 
the (messy) statistics at distances k <^ 6r <^ dt. We want to plot u = R/dt versus N/gt 
where R denotes either the diameter (full symbols) or the radius of gyration (empty 
symbols). We still need to fix the persistence length dependence of df. For simplicity, we 
suppose that dt is only weakly affected by a and set arbitrarily dt = I, i.e., we neglect 
any stiffness dependence of dt itself. Hence, gt oc l/gk and v = Ng^. This assumes 
that the effect of a is to increase R/dt (allowing the asymptotic regime to be accessed 
for smaller A^) by increasing R at (nearly) constant dt. This idea is very much the 
opposite of what we tried in Fig. ^ where we had v = N/gk and stiffer chains with fewer 
statistical units were considered as being effectively "shorter" . As shown in Fig. ^ this 
simple new proposal is successful although a weak dependence of dt on a cannot be ruled 
out. Moreover, it is self-consistent: The local statistics is indeed well described by the 
Gaussian behaviour z/i ^ 1/2. This might be interpreted as the a ^ limit of the 
CD proposal. The data collapse over two orders of magnitude in v justifies a posteriori 
the neglect of a stiffness dependence of dt. At large distances data points for different 
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persistence lengths fall together on one slope of 1^2 = 1/3 which is consistent with the LA 
picture. Note the later crossover of the radius of gyration onto the z/2-slope. This might 
be related to the fact that Rg probes smaller distances. 

The above proposal may be in principle generalised to incorporate excluded volume 
effects on intermediate scales A; ^ (5r ^ ^ where the statistics is governed by the exponent 
z/q = 0.59 for dilute rings. Note that the strong dependence of ^ on the persistence length 
is checked by the small difference between uq and z/i. (We plan to consider this problem 
in a subsequent study on density effects [Q.) This offers an effective simplification for 
the molten chains (0 = 0.5) of interest in this study: Both length scales ^ and k are of 
same order and we are not able to separate the effects anyway as we are going to show 



now. 



F. Local conformational properties: Structure factor 

So far we have considered mainly global properties like the ring diameter R^ or the 
ring area. We wish now to characterise the local conformational properties by measuring 
the single chain structure factor S{q) and the mean-square length (Rl), both introduced 



in Sec. |T|, and the differential fractal dimensions df{q) associated with both quantities. 
We will first consider the fractal dimension obtained from S{q) as presented in Fig. ^ and 
then compare this in Sec. [IV F| below) with df obtained from (R^) and shown in Fig. |10 



The differential fractal dimension for our longest (A^ = 1024) rings in the melt (0 = 
0.5) is depicted in Fig. |^(a) and compared with their dilute counterparts. The Bragg 
peak at 6r ^ I is much more pronounced than for the dilute chains (not fully shown). 
We note again small, but distinct cr-dependent packing effects. Not surprisingly, at very 
large distances the structure is again well described by the Guinier expansion (see section 
pnp and df{q) vanishes smoothly. 

Obviously one wants to understand the behaviour between the (featureless and trivial) 
Guinier and the (non- universal) Bragg part. To stress this regime we have plotted in 
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Fig. ^(b) the differential fractal dimension df{q) versus Sr = 271 /q chopping off distances 
5r < I and the Guinier part. Data are for flexible (cr = 0) and semi-flexible {a = 3) 
rings and linear chains in the melt. For small distances where df < 1 the data for all 
df are very similar. (Note the small a-dependent packing effects mentioned above.) At 
df{6r ^ A) = 1 the semi-flexible systems branch off the much steeper line for flexible 
chains (top line). This applies also to the data for cr = 1 and a = 2 (not shown). We 
stress that the difference between different a is gradual and that systems of flexible rings 
do not constitute a singular limit. Because of the small blob size (i.e., the very strong 
interactions with the surrounding rings) we are unable to separate rigidity and excluded 
volume effects: Semi-flexible rings do not show the "shoulder" at df ^ 1 we saw for dilute 
rings (dotted line in Fig. H(a)). 

At these small distances, smaller than the blob size ^ [|15[, linear chains behave exactly 



like their ring counterparts of the same flexibility. This proves, as assumed above, that 
on small scales the a-effects are the same for both architectures and justifies the use of 
the persistence length b (and the associated values Qk and k) obtained from linear chains. 
It also shows that it was natural to attempt in ref. |Q a finite-size scheme in terms of 
the correlation hole, at least when one assumes the CD picture to hold, i.e. R being the 
only length scale. 

However, the rescaling in Fig. ^ did not work for molten rings (in contrast to linear 
chains). This is due to effects caused by the ring closure and topology which intervene 
at larger distances 6r ^ ^ (and in particularly at 6r ^ df) : The fractal dimension df for 



rings continues to rise while the linear chains become roughly Gaussian, i.e. df ^ 2 [p4 
The more flexible the rings, the more this increase becomes suppressed (see the flexible 
chains of mass A^ > 256 in Fig. |^(b)). This gives the more rigid systems a chance to catch 
up (in df{q)) the more flexible ones — and this after having consumed fewer monomers 
on short scales. 



The above discussion is only valid if the chains are actually large enough i.e. 6r ^ R 



9 
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or S{q) ^ A^, to fall on the asymptotic (cr-dependent) curves. The closure constraint 
(rather than the topology) forces smaller rings (e.g., iV = 64 in Fig. §(b))), to form rather 
dense globules of blobs. 

For (asymptotically) long rings there are now two possible scenarios: Either all the a- 
lines merge again at large enough distances or they become parallel. (We reject crossing 
of the df{a)-\mes as unphysical.) Unfortunately we are at present unable to decide 
unambiguously between these two. We believe it however more likely that all the lines 
eventually merge and that the asymptotic properties of ring polymers are independent 
of the persistence length b p5|. Hence, for large distances 6r > 6r*{a) the evolution of 



df{6r) becomes again independent of a. It is tempting in view of the scaling in Fig. ^ to 



set Sr* = dt\2t 



G. Local conformational properties: Contour distance i?„ 

We consider finally the average distance i?„ = {{Ri — Ri+nYY^'^ between the ends of 
a ring contour segment of length n. Similarly to the discussion above for S{q) we define 
a differential fractal dimension l/df = d\og{Rn) / d\og{n) . Note that S{q) and (R^) are 
not just simply Fourier transforms to each other. They contain different information: the 
S{q = 27r/i?„) monomers of the reference ring within a volume of radius ~ i?„ around an 
arbitrary monomer of the same ring comprise also monomers far away the ring contour, 
i.e. much larger than n. This contribution is certainly small for linear chains (hence 



S{q) (X n 1T2|), but gets more and more important while the rings become more compact 
(with increasing a and A^). 

We have plotted Rn in the main part of Fig. |10| versus n < N/2 for various ring 
lengths mass as indicated in the figure. Both flexible (cr = 0) and semi-fiexible {a = 3) 
are included. The lines of the stiffer rings are stronger curved. For n <^ N/2 all curves 
collapse on a chain length independent (but a dependent) master curve. Obviously, 
in the limit of n — ;> N/2 the ring closure forces Rn to level off towards Re, i.e. l/dj 
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has to vanish. The inset shows the differential fractal dimension df for our longest 
rings A^ = 1024 for three different persistence lengths. Indeed, df for small n increases 
with 0". For short flexible rings the fractal dimension takes more or less the Gaussian 
value. Consistent with Fig. |^ we are unable to separate clearly the different length scales 
and observe broad crossover lines. As mentioned above, df has trivially to diverge for 
n ^ N/2 and data points n > N/A are influenced by this upper cut-off. However, we 
again see that stiffer chains show more compact scaling than flexible ones. Qualitatively, 
the a = 3-curve is even consistent with the df = 4-window predicted for Gaussian LA 

i- 

In brief, the advantage to use stiff chains is due to an additional length scale dt oc A^''. 

Flexible chains, very compact on short length scales (see the snapshot on the left of 

Fig. 0), "waste" there a large number gt of monomers. Systems with N < gt (that is, 

Rg <^ dt) feature only one characteristic size, their own size, and show CD scaling (with 

a = 0). Semi-flexible chains, however, using less monomers on short distances below dt 

can explore larger distances and become more compact in the scaling sense: They do this 

by having a lower ceiling on the local density p of a single ring, thereby enhancing the 

topological interaction with their neighbours (thus increasing df). 

V. DYNAMICS OF RINGS IN THE MELT 

The differences in the ring statistics should strongly influence their dynamic be- 
haviour. But as pointed out in ref. 0], there is some experimental evidence that the 
dynamics of rings are similar to their linear counterparts, at least up to the largest 
molecular masses that can readily be obtained. Our simulation data, presented below, 
confirms this. 

We characterise the dynamics by measuring three different mean-square displace- 
ment functions describing the motion of monomers gi{t) = ([R„(t) — R„(0)]^) in the 
laboratory frame, the motion of monomers in the center-of-mass frame of a given ring 
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g2{t) = {[^n{t) — Rc.m.(t) — R„(0) + Rc.m.(O)]^), and the motion of the center of mass 
itself (y'3(t) = ([Rc.m.(^) ~ Rc.m.(O)]^)- ^^ shown in Fig.|Tl|, the mean square displacements 
for flexible rings collapse onto chain length independent master curves when the mean 
square displacements are rescaled by the radius of gyration (-Rp, and the time by the 
characteristic relaxation time {B?) / D^- Here D^q is the diffusion coefficient (presented 
in Fig. |l^) obtained from the asymptotic behaviour of the center of mass motion g^{t)/t 
(bold line of slope 1 on the right). 

Closer inspection shows, however, that the collapse is not perfect. This is in agreement 
to what was observed by Brown and Szamel P]. In any case, this scaling is not born 



out for the semi-flexible (a = 3) systems of mass A^ = 128 depicted in Fig. |TT|. This 



does not come as a surprise in view of what we have described above in section |V 
about the scaling of conformational properties. It is worth noting than in inspection of 
the differential slope of the motion of monomers gi (t) (not shown) shows (a small, but 
distinct) region with slope 0.32 instead of the classical Rouse-like anomalous diffusion 
exponent ?« 1/2. 

In Fig.|T2| we present the self-diffusion constant of linear chains (empty symbols) and 
rings as a function of their size A^. The circles correspond to isolated swollen chains. 
The diffusion constant is the same for linear chains and rings. It scales like Dn ~ 1/A^. 
Our Monte Carlo simulations do not of course incorporate hydrodynamic interactions 



and the dynamics is expected to be Rouse-like |jT3|. The squares correspond to flexible 
chains in the melt = 0.5. It is important to note that rings are always /aster than their 
linear counterpart of same mass. As shown in the flgure, the diffusion coefficient are well 
fitted by (effective) power laws: Dn oc A^~^-^^ for rings and D^ oc A^~^'^ for linear chains. 
As was emphasised in ref. Q] both linear chains and rings scale with the size of their 
correlation hole Dat oc R^^. The diamonds in Fig.|l^ display the results of semi- flexible 
rings {a = 3) at melt density (p = 0.5. The dependence of the diffusion constant for 
the larger ring sizes studied is stronger than for flexible chains and obeys the apparent 



20 



relation D^ ~ A^~^'^^. This is roughly consistent with the data of Brown and Szamel 
who found D ~ A^~^-^^ for large flexible rings 0. Both sets of simulation data are likely 
to be affected by crossover effects. 

This observation is in accord with theoretical calculations, which calculate the diffu- 
sion constant of rings via the motion of kinks along the contour. This diffusion of kinks 
along the molecules is also the dominant relaxation mechanism in the reptation of linear 
chains. In this sense, the dynamics of (more and more entangled) rings is the analog of 
reptation in melts of linear molecules. We do not find any signature of arm retraction in 
the motion of ring polymers. The latter would result in an exponential decrease of the 
diffusion with ring size like in star polymers ||27|| . 

In the (suspected) LA limit, the ring conformations should possess an hierarchical 
structure. Studying the motion of a single ring in a network of (explicitly) fixed obstacles, 
Obukov et al. suggested that the outer arms of a LA can rearrange much faster than the 
inner structure 0]. This results in a broad range of relaxation times. We investigated this 
by monitoring the distribution of mean square displacements of monomers for flexible 
rings {N = 256 and a = 0) and semi-flexible ones (A^ = 256 and a = 3) after a time 
interval At ^ 2.3 10^ MCS. This corresponds to a time scale which is shorter than the 
diffusion time D^/ Rg of the ring. Note that both curves correspond to roughly the same 
time period and the same mean value of the monomer displacement. The center of mass 
diffusion constant differs, however, by an order of magnitude. Hence, the major difference 
in the dynamic behaviour does not stem from a difference in the local dynamics (e.g. a 
slight decrease of the acceptance ratio of the Monte Carlo moves due to the additional 
bond angle potential), but reflects the interplay between conformations/topology and 
the dynamics. Notably, the monomer mean square displacements for flexible rings are 
Gaussian distributed, while the distribution of semi-flexible systems is non-Gaussian. 
This is consistent with simulations of rings in a network of fixed obstacles [0]. 

From these different probes of the ring dynamics we conclude that the dynamics for 
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flexible and semi-flexible chains are qualitatively different, at least for the chain lengths 
we were able to probe. There appears to be clear evidence for the occurrence of a second 
length and time scale for our semi-flexible chains (which have R > dt), but not for our 
flexible ones (which have R < dt). Qualitatively we conclude again that stiffer rings are 
closer to the LA limit than their flexible counterparts. However, even for the largest ring 
size and stiffness our data are still affected by the small number of arms and we expect 
pronounced flnite size effects to the asymptotic behaviour for large rings in the melt. An 
investigation of the dynamics of even larger rings {N > 256) is unfortunately beyond our 
computational facilities at present. 

VI. CONCLUSIONS 

In summary, we have presented extensive Monte Carlo simulations of semi-dilute 
solutions and melts of flexible and semi-flexible ring polymers. The rings are neither 
knotted with themselves nor concatenated with each other. 

In order to have a reference with systems where the non-linkage constraint does not 
play a role we have briefly considered dilute solutions of rings. They appear to be ex- 
tremely similar to their linear counterparts, both with respect to statical and dynamical 
properties. The (remaining) topological unknottedness constraint seems not to be perti- 



nent |2T|, at least not for the chains sizes (A^ = 1024) we have considered here. The ring 
extension scales like R ~ N'^ with z/ ^ 0.59. This was also made evident from the study 
of the differential fractal dimension df{q) obtained from the static structure factor S{q). 
Dilute rings do not differ from linear chains with regard to the influence of flnite persis- 
tence length (at least for large enough rings). The dynamics is Rouse-like Dn ~ A^~^, 
because our Monte Carlo simulations ignore hydrodynamical interactions ( "free draining 
limit"). In short, from the practical point of view topology is irrelevant. 

This is dramatically different in the high density limit (0 = 0.5) with strongly overlap- 
ping entangled rings where topological constraints tend to squeeze the rings into relatively 
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compact objects. We have varied the stiffness of the rings so as to tune the overlap be- 
tween different rings. Due to the rather compact structure of the molecules in the melt, 
increasing the stiffness is much more efficient than increasing the ring size. Essentially, 
stiffer chains "waste" less monomers on short distances and have more monomers left to 
meander through the topological constraints imposed by neighbouring rings. Indeed the 
effective Flory exponent i^iN) obtained from the high chain length behaviour of chain di- 
ameter and radius of gyration shows a strong effect with regard to the persistence length, 
decreasing from u ^ 0.4 for flexible chains (cr = 0) to z/ ~ 1/3 for our stiffest systems 
(cr = 3). Chain stiffness allows a reduction in u to more "compact" values by increasing 
(at a given u) the overlap parameter p oc b^N^'^^^. 

Rings with topological constraints do not follow the classical one parameter scaling 
with 0/0* for linear chains where the size of the (dilute) chain of mass A^ and stiffness 



a sets the only relevant length scale |T^. This is in disagreement to the fundamental 
assumption of the CD picture 0. In order to scale the chain length R = R{N, a) we were 
forced to assume an additional chain length independent length scale df oc gt'^^ . Supposing 
a weak cr-dependence of dt and choosing ui self-consistently this yields a satisfactory 
data collapse and ui ^ 1/2 as demonstrated in Fig. |^. For small chains, R <^ df, 
the topological interactions are weak. The rings resemble their linear counterparts and 
behave effectively like closed Gaussian chains of blobs. This regime is consistent with the 
a — i> limit of the CD scenario, i.e. due to the non-linkage constraint the free energy 
of a reference ring increases only like log(p). Rings larger than dt appear to be governed 
by the topological interactions and are more compact (in the scaling sense of smaller 
u). They are well characterised by an effective Flory exponent z/2 = 1/3. This scaling 
scenario is broadly consistent with the concept of lattice animals (LA) within a network 
of topological obstacles created in a self-consistent manner by surrounding rings [0,0. 
These LA appear to be made of locally Gaussian chain pairs at short distances. 

Our (more ambitious) discussion of the differential fractal dimension dj (obtained 
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from the static structure factor S{q) and/or the spatial distance Rn{n) between a contour 
segment of length n) reveals very rich and broad crossover effects: excluded volume 
^ versus persistence length fc, excluded volume versus ring closure, ring size R versus 



topological length scale dt- As shown in Fig. y and Fig. |T0|, it is not possible to separate 
the different length scales unambiguously and to disentangle their physics [^. In view 
of the (restricted) range of parameters (A^ < 1024, cr < 3,0 = 0.5) we are able to 
simulate, this does not come as a surprise. Much more surprising is the success and the 
simplicity of the scaling scenario of Fig. ^ for the global chain R size described in the 
paragraph above. There, all the intricate short range physics was cast in one effective 
exponent vi for all chain and persistence lengths. It just turns out that the Gaussian value 
Vi = 1/2 (i.e. a — > 0) fits particular well the data. Similarly, while 1/2 = 1/3 is certainly 
the asymptotic value, this does not exclude the possiblity of an intermediate window 
with z/2 = 1/4. Indeed the differential fractal dimensions for our largest and stiffest 
configurations clearly exceed df = 3. This is in favour of an ultra-compact transient 
which should then eventually become evident also in the ring sizes for even larger chain 
lengths A^ than we are at present able to simulate. 

Additional evidence for the crossover to a strongly entangled regime characterised by 
an additional length scale dt comes from our brief investigation of the ring dynamics. 
The scaled time dependence of the monomer displacements differs from the master curve 
in the non-entangled regime, and the monomer displacements at times smaller than 
the relaxation time are non-Gaussian distributed as expected for LA |0|. The diffusion 
constant for chains in the melt scales like Dj^i ~ A^"^-^^ for flexible chains. This decrease 
to Dn ~ A^-i-68 for semi-flexible systems (a = 3). Again these exponents are presumably 
only effective values due to a broad crossover between unentangled and entangled regime. 
The similarity of the dynamics of melts of rings and linear chains suggests that our 
observations for rings might also be pertinent to the dynamics of linear chains. Indeed, it 
is tempting to relate the topological constraints which lead to the lattice animal behavior 
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for rings to the entanglements in linear chains. While the topological interactions do not 
influence the static conformations of linear chains, however, rings offer the additional 
possibility to investigate the effect of topology in the static behavior. 

In any case, our simulation data cover only the onset of LA behaviour and our es- 
timates for the scaling functions of the conformational statistics and the dynamics are 
likely to be subjected to corrections due to the very small number of arms. In future 
we plan to corroborate further the discussion of the dynamical properties (increasing 
the number of statistical segments) and to investigate the static and dynamical scaling 



properties with regard to the monomer density |14 
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TABLES 
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-0.668 


2.703 


4.24 
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TABLE I. Persistence length dependent properties for dilute {(j) = 0) rings and linear 
chains: Mean cosine of the bond-bond angle (cos(0)), mean bond length /, effective bond 
length of linear chain 6g and rings 6g , acceptance rate A and the diffusion constant ND^. 
These values characterise the asymptotic behaviour of long linear chains as well as rings. 
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0.1136 


0.36 
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-0.676 


2.593 


5.3 


0.0983 


0.33 



TABLE II. Persistence length dependent properties for rings and linear chains at high 
volume fraction (p = 0.5: Mean cosine of the bond-bond angle (cos(B)), mean bond length I, 
effective bond length b^ for linear chains, acceptance rate A and the differential Flory exponent 
Ueff fitted over the four largest chains available (see Fig. y). 
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TABLE III. Ring radius of gyration Rg, diameter i?e and (reduced) diffusion constant 
A^DAflO^ for different ring sizes N and stiffnesses cr at </> = 0.5. 
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FIGURES 




(c) 




FIG. 1. Sketch of topological constraints and their effects, (a) The top chains present 
allowed conformations of unknotted and (more importantly) unconcatenated rings. The rings 
cannot turn into knotted configurations (not shown) or the concatenated configurations on 
the bottom. It is this non-linkage constraint which tends to squeeze a ring in semi-dilute 
solutions and melts. The interplay between the topological constraint and the effective bond 
length b (in contrast to the monomer size a) is studied here, (b) One ring in a network of 
fixed topological constraints or obstacles (squares) imposing a strongly entangled and compact 
conformation, (c) An equivalent Lattice Animal (LA): In the high overlap limit rings are 
expected to behave effectively like LA. This is mean-field picture for a reference ring within a 
self-consistent topological network imposed by the surrounding rings. 
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FIG. 2. Configurational snapshots of rings of mass N = 256: dilute flexible ring in the 
middle, flexible (a = 0) ring in the melt {(p = 0-5) on the left and corresponding semi-flexible 
{a = 3) coil on the right. The overlap has increased by factor three between the flexible and 
the semi-flexible chain in the melt. To obtain a similar effect by increasing the chain mass we 
would have to increase A^ by at least two orders of magnitude. 
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FIG. 3. Diameter i?e versus chain length N for different a for dilute (open symbols) and 
dense (full symbols) rings. The dilute rings are characterized by the same Flory exponent 
V = uq ~ 0.59 as their linear counterparts (not shown). For large enough chains this is 
independent of a. For the dense chains the effective fractal dimension df = 1/u becomes 
much larger, depending strongly on the chain rigidity: flexible chains (a = 0) compare well 
with V = 0.4 (dashed line), semi-flexible chains {a = 3) with v = 1/3 (solid line). 
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FIG. 4. Reduced diameter u = Re/k versus reduced chain mass v = N/gj^ for dilute {(p = 0, 
open symbols) and dense {(p = 0-5, full symbols) linear chains (top two lines) and rings (bottom). 
The chain flexibility a changes as indicated in the figure. The lines indicate (effective) exponents 
(from top to bottom): v = 0.59 for dilute linear chains, v = 0.5 for dense linear chains, u = 0.59 
for dilute rings, i' = 0.4 for dense, flexible (a = 0) rings and z^ = 1/3 for semi-flexible {a = 3) 
rings in the melt. Note that the dense rings do not collapse on one master curve as their linear 
counterparts. 
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FIG. 5. Bond-Bond correlation function < eoe„ > for dilute (</> = 0) and molten {(p = 0.5) 
rings for o" = (connected by thin lines) and o" = 3 (connected by fat lines) for various chain 
masses as indicated in the figure. At high chain overlap (high density and rigidity) we find a 
pronounced negative dip in the correlation function. 
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FIG. 6. Differential fractal dimension df{q) = —dlog{S (q)) / dlog{q) versus Sr = 2iT/q. (a) 
Dilute and molten rings of mass N = 1024. The line indicates an ideal Gaussian ring of same 
radius of gyration as a flexible and dilute ring, (b) Rings in their melt {(j) = 0.5) for N = 1024 
(circles), N = 512 (squares), N = 256 (diamonds) and A^ = 64 (triangle). The dotted line is 
a comparison with linear chains {N = 256, (f> = 0.5). Within the first blob {6r < ^(cr)) linear 
chains and rings behave identically. 
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FIG. 7. Single chain density distribution with respect to its center of mass for dense systems 
(j) = 0.5 with flexible rings (cr = 0) versus reduced distance r/Rg. The density at the center of 
the coil decreases {Rg /N oc N^^'"^ in agreement with the effective exponent from Fig. ^, but 
remains always larger than 0.4(/>. The inset presents the same for semi-flexible rings [a = 3); 
the (rescaled) distribution density becomes chain length independent for A^ < 64 (density at 
the center « 0.25(/)). 
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FIG. 8. Reduced area {\a\) / {B^) versus N for dilute rings (open symbols) and molten rings 
(full symbols) for different a. Note the logarithmic, but systematic correction to the scaling for 
the latter. 
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FIG. 9. Two scaling attempts for diameter (full symbols) and radius of gyration (empty 
symbols) at different flexibilities as indicated in the figure. In the inset we try the traditional 
scaling u = R/Rq versus v = (p/cp* ~ R^cp/N. The poor scaling evidences additional length 
scales. In the main figure we explore the possibility of a length scale dt assumed to be inde- 
pendent of both mass N and persistence length a. We plot u = R/dt versus v = N/gt where 
we set arbitrarily dj = 1 and fix gt/gk = (dt/k)'^ self-consistently with the observed effective 
Gaussian behavior at short distances. All data points collapse. The slopes indicate Gaussian 
behavior for R < dt (i^i = 0.5: solid line) for R < dt and compact LA (z/2 = 1/3: dashed line) 
for larger distances. 
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FIG. 10. Distance Rn between a chain contour of length n for <j) = 0.5. Main figure: Flexible 
and semi-flexible {a = 3) chains of various chain length N as indicated in the figure. Inset: 
Differential fractal dimension df{n) = \/u[n) for N = 1024 and cr = (dotted line), a = 1 
(dashed line) and a = 3 (solid line). 
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FIG. 11. Reduced mean-square displacements versus reduced time for two rigidities and 
different chain length as indicated in the figure. The scaling works reasonably well for flexible 
chains (a = 0). However it is inconsistent with the semi- flexible (a = 3) configurations. 
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FIG. 12. Diffusion constants ND^ reduced by tlie acceptance rate A versus N for rings 
(full symbols) and linear chains (empty symbols). We consider dilute, flexible (o" = 0) chains 
(circles) and flexible (a = 0, squares) and semi-flexible {a = 3, diamonds) chains in the melt 
((/> = 0.5). The lines correspond to effective power laws. 
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FIG. 13. Probability distribution of mean-square displacements of monomers at short 
times. Flexible rings (0 = 0.5, cr = 0): Gaussian distribution of monomer displacements at 
QDAt/R^ = 0.42. Semi-flexible rings {4> = 0.5, a = 3): deviations from the Gaussian behavior 
at 6DAt/R^ = 0.02 indicate abundance of slow monomers (presumably) within the "trunk" of 
the lattice animal tree. The solid line represents the Gaussian distribution. 
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